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In order to describe the exotic nuclear structure in unstable odd-A or odd-odd nuclei, the de- 
formed relativistic Hartree Bogoliubov theory in continuum has been extended to incorporate the 
blocking effect due to the odd nucleon. For a microscopic and self-consistent description of pairing 
correlations, continuum, deformation, blocking effects, and the extended spatial density distribution 
in exotic nuclei, the deformed relativistic Hartree Bogoliubov equations are solved in a Woods-Saxon 
basis in which the radial wave functions have a proper asymptotic behavior at large r. The for- 
malism and numerical details are provided. The code is checked by comparing the results with 
those of spherical relativistic continuum Hartree Bogoliubov theory in the nucleus '°O. The prolate 
deformed nucleus ‘°C is studied by examining the neutron levels and density distributions. 


PACS numbers: 21.60.-n, 21.10.-k, 21.60.Jz 


As one of the most fascinating phenomena found in 
exotic nuclei, the nuclear halo contains many new and 
interesting features, such as a largely extended spatial 
density distribution, dilute and pure neutron matter, the 
Borromean property, and the coupling between bound 
states and the continuum [1]. In the past, several rela- 
tivistic Hartree Bogoliubov theories with or without the 
Fock term have been developed for a self-consistent de- 
scription of spherical halo nuclei [2-6]. In order to in- 
vestigate halo phenomena in deformed nuclei in a mi- 
croscopic and self-consistent way, a deformed relativistic 
Hartree Bogoliubov (RHB) theory in continuum has been 
developed for even-even nuclei [7-9]. Within this model 
a decoupling in shape between the core and the halo has 
been predicted in some deformed nuclei close to the neu- 
tron drip line, e.g, in 4744Mg [8, 9]. 

In order to describe the exotic nuclear structure in un- 
stable odd-A or odd-odd nuclei, the blocking effect of one 
or several nucleons has to be taken into account. In the 
present Letter, the deformed relativistic Hartree Bogoli- 
ubov theory in continuum [7-9] is extended to incorpo- 
rate the blocking effect due to an odd nucleon. In such a 
way, pairing correlations, continuum, deformation, block- 
ing effects, and extended spatial density distributions in 
exotic, odd-A or odd-odd nuclei can be taken into ac- 
count microscopically and self-consistently. 

In order to treat pairing correlations, the quasi-particle 
concept is adopted and the ground state of an even-even 
nucleus |®) is represented as a vacuum with respect to 


quasi-particles [10], 
Gx|®o) =0, forall k=1,...N, (1) 
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where N is the dimension of the quasi-particle space and 
the quasi-particle operators pl, By are obtained by the 
Bogoliubov transformation from the particle operators 


d,c, 
B c U V* 
(aar) W= (oy) (2) 


The matrix W of the coefficients U and V is unitary 
which guarantees that the quasi-particle operators satisfy 
the same anti-communication rules as the particle opera- 
tors [10]. Starting from the bare vacuum |0}, the ground 
state of a system with an even number of particles |®o) 
can be constructed as, 


|®o) = | [ 210), (3) 
k 


where k = 1,...N. For an odd system, in practice, the 
corresponding ground state can be constructed as one 
quasi-particle state |®,), 


|®1) = 6] |o) = Bf [ [ 2410), (4) 
k 


where pi corresponds to the quasi-particle state with the 
lowest quasi-particle energy. In other words, the one 
quasi-particle state |1) is the vacuum with respect to 
the set of quasi-particle operators (681, ..., 8) with 


CF). BHPian fyi (5) 


and the exchange of the operators gt +> 6, forms a new 
set of quasi-particle operators (6/,...,8'y, i',..., BA), 
which corresponds to the exchange of the columns 
(Un, Va) <> (Vi, UÑ) in the matrix W. That is, the 
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blocking effect in the odd system can be realized by ex- 
changing the creator fal with the corresponding annihila- 
tor 6; in the quasi-particle space. Accordingly, the block- 
ing effect in a multi-quasi-particle configurations can be 
treated. 

The covariant density functional theory has been ap- 
plied to describe successfully nuclear structure over the 
entire periodic table [11-15]. The relativistic Hartree Bo- 
goliubov equations [13-15] for the nucleons read, 


jane Pee (Ur) =E a (6) 


where Ep is the quasiparticle energy, À is the chemical 
potential which guarantees the proper average particle 
number and Ap is the Dirac Hamiltonian, 


hp(r,r’) =a-p+Vi(r)+8(M + S(r)). (7) 
The scalar and vector potentials 
S(r) = goo(r), (8) 


= Ar), (9) 
depend on the scalar field ø as well as the time-like com- 
ponents w, p°, and A? of the iso-scalar vector field w, the 
3-component of iso-vector vector field p and the photon 
field A. 

For a fully paired and axially symmetric deformed sys- 
tem with the time reversal symmetry, the projection of 
the total angular momentum on the symmetry axis Q is 
conserved and each single particle state has a degeneracy 
of two. The RHB equation (7) can be reduced to half 
dimension M = N/2 and can be decomposed into de- 
generate blocks with quantum numbers +2 or —Q. The 
corresponding density and abnormal density matrix have 
dimension M and read 


V(r) = gow? (r) + 9pT3P(r) +e 


Pmxm — V*VT, (10) 
k = V*UT, (11) 


MxM 
where V and U are the coefficients in the corresponding 
subspace. 

For an odd system with the ky-th level blocked in the 
+Q subspace, the time reversal symmetry is violated and 
there appear currents in the system. These currents show 
an axial symmetry, i.e., Q remains a good quantum num- 
ber, but the quasi-particle energies are no longer degener- 
ate for the two subspaces, because the subspace with +9 
contains the odd particle and the corresponding subspace 
with —Q contains an empty level. Therefore, in princi- 
ple, we have to diagonalize the RHB equation twice, one 
for the subspace with +Q and the other for the subspace 
with —Q. Using the equal filling approximation which is 
usually made [16, 17], we neglect currents and average 
in a statistical manner over the two configurations of a 
particle in the +Q space and a particle in the —Q space. 
The corresponding currents cancel each other and in this 
way we obtain in each step of the iteration fields with 


TABLE I. Ground state properties of 19O from deformed 
RHB and RCHB calculations. Rn, Rp and Rs refer to neutron, 
proton and total root mean square radii in unit of “fm”. Epart, 
Es, BS", Ea, He", Ep, Ecoul, Epair and Erot refer to energies 
of particle, o field, non-linear term of o field, w field, non- 
linear term of w field, p field, Coulomb field, paring energy 
and total binding energy in unit of “MeV”. 


deformed RHB RCHB 
Ry 2.83 2.83 
Rp 2.57 2.57 
Ry 2.72 2:72 
Epart —445.60 —445.73 
Es 2259.26 2260.18 
Be" —77.14  —77.18 
Bus —1880.10 —1880.83 
Epo 32.76 32.77 
Ep —3.15 —3.16 
Ecos —17.10 —17.10 
Epair —3.82 —3.82 
Erot —144.80 —144.78 
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FIG. 1. (Color online) Neutron single particle levels in the 
canonical basis and their occupation probability v? for O 
from the deformed RHB and RCHB calculations. 


the time reversal symmetry. In practice we average the 
density matrix p and symmetrize the abnormal density 
k in two subspaces and replace the two densities by 


1 * * 

p = Puxm + 5 (Ur, UT E Vig, Vin) ) (12) 
1 * * 

K = BmxmMm ` 5 (Un, Ve" at Vi. Ug, ) ’ (13) 


where Vk, and Uk, are column vectors in the matri- 
ces V and U corresponding to the blocked level. Note 
that Eqs. (12) and (13) are the same as those given in 
Refs. [16, 17] if one considers the time reversal symmetry 
and that the dimension of the densities here is M = N/2. 

In the following applications we use the density func- 
tional PK1 [18] and a zero range density dependent pair- 
ing force [4, 9]. The deformed RHB equations are solved 
in a spherical Dirac Woods-Saxon basis [19] and the de- 
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FIG. 2. (Color online) Neutron density profiles of ‘*C and °C 
in deformed RHB calculations with the parameter set PK1. 


formed potentials and densities are expanded in terms 
of the Legendre polynomials P\(cos@) (for details see 
Refs. [7—-9]). In order to check the accuracy of the de- 
formed code with the blocking effect we investigate the 
spherical nucleus 190 and compare the results of the de- 
formed RHB code allowing only the spherical compo- 
nents of the fields, i.e., A = 0, with those obtained in 
the spherical relativistic continuum Hartree Bogoliubov 
(RCHB) theory [4]. The ground state properties from 
the deformed RHB and the RCHB calculations are given 
in Table I. It is clearly seen that both calculations agree 
well with each other. For example, the difference between 
the total binding energies is about 0.02 MeV, which cor- 
responds to an accuracy of about 0.01%. 

In order to check the details of the deformed RHB cal- 
culation, we show in Fig. 1 the neutron single particle 
levels in the canonical basis for 19O in comparison with 
the spherical RCHB results. The length of each level is 
proportional to the occupation probability. In the spheri- 
cal RCHB calculation, the blocked orbital is 1d5/2. In the 
deformed RHB code, a spherical solution is enforced by 
allowing only A = 0 component in the Legendre expan- 
sion mentioned earlier and therefore the three sublevels 
of the 1d5/2 orbital with Q7 = 1/2*, 3/27, and 5/2* are 
degenerate. The present results in Fig. 1 are obtained 
with the level Q7 = 1/2* blocked in the deformed RHB 
code, but blocking any of these three levels gives the same 
results. From Fig. 1 one finds a good agreement between 
the results from the deformed RHB and RCHB calcula- 
tions. 

After checking the code in the spherical case we now 
turn to an application of the deformed RHB theory for 
an exotic and deformed nucleus with an odd number of 
particles. The deformed halo candidate nucleus °C [20] 
is studied with the parameter set PK1 [18]. In contrary to 
the spherical even-even nucleus 4C, the deformed RHB 
theory predicts a deformation 8 = 0.25 for the odd-A 
nucleus 15C. The neutron root mean square radii of 14C 
and 15C are 2.56 fm and 2.79 fm respectively. 

The deformed RHB calculations lead for 15C to a de- 
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FIG. 3. (Color online) Neutron single particle levels in the 
canonical basis and their occupation probability v? for °C in 
the deformed RHB calculation. 


formed density distributions with an axial symmetry. In 
Fig. 2 the neutron densities along the symmetry axis 
Pn(rz,r1 = 0) and perpendicular to the symmetry axis 
Pn(rz = 0,r1) are plotted as dashed and dashed-dotted 
lines. We also include the spherical density for 14C as 
a reference by a solid line. It is interesting to see that, 
in the direction perpendicular to the symmetry axis, the 
neutron densities are almost the same for 14C and °C 
at least when r} < 6 fm. Along the symmetry axis the 
neutron density of 15C extends much further than that 
of 14C. This is partly due to that 15C is prolate and the 
weakly bound 1/2* level is occupied. 

The single particle levels of C in the canonical basis 
are plotted in Fig. 3. Since it is a deformed nucleus, the 
1p3/2 orbit is split into two levels with Q7 = 1/27 and 
Q7 = 3/27 respectively. There is one neutron occupying 
in Q7 = 1/27 level near the threshold. The occupation 
probability v? = 0.5, indicating that it is averaged over 
the two configurations with Q7 = +1/2*. Because of 
the deformation, this level is a mixture of the spherical 
orbits 1d5/2 (62%) and 28/2 (36%). The weakly-bound 
feature and the relatively large s-wave component of this 
level results in that the neutron density of °C extends 
further along the symmetry axis. 

In summary, the blocking effect due to an odd nucleon 
is incorporated in a deformed relativistic Hartree Bogoli- 
ubov (RHB) theory in continuum in order to describe 
odd-A or odd-odd exotic nuclei. The formalism is briefly 
presented and the numerical checks are carried out by 
comparing the results of the deformed RHB code for the 
spherical nucleus !°O with those obtained from the spher- 
ical RCHB code. As a first application, the deformed 
nucleus !°C is studied by examining the neutron levels, 
deformation as well as density distributions along and 
perpendicular to the symmetry axis. Along the symme- 
try axis the neutron density of !°C extends much further 
because !°C is prolate and the weakly bound valence level 
with Q7 = 1/2* has a relatively large s-wave component. 
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